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We present an analytical theory for the gate electrostatics and the classical and quantum capaci- 
tance of the graphene nanoribbons (GNRs) and compare it with the exact self-consistent numerical 
calculations based on the tight-binding p-orbital Hamiltonian within the Hartree approximation. 
We demonstrate that the analytical theory is in a good qualitative (and in some aspects quanti- 
tative) agreement with the exact calculations. There are however some important discrepancies. 
In order to understand the origin of these discrepancies we investigate the self-consistent electronic 
structure and charge density distribution in the nanoribbons and relate the above discrepancy to the 
inability of the simple electrostatic model to capture the classical gate electrostatics of the GNRs. 
In turn, the failure of the classical electrostatics is traced to the quantum mechanical effects leading 
to the significant modification of the self-consistent charge distribution in comparison to the non- 
interacting electron description. The role of electron-electron interaction in the electronic structure 
and the capacitance of the GNRs is discussed. Our exact numerical calculations show that the 
density distribution and the potential profile in the GNRs are qualitatively different from those in 
conventional split-gate quantum wires; at the same time, the electron distribution and the poten- 
tial profile in the GNRs show qualitatively similar features to those in the cleaved-edge overgrown 
quantum wires. Finally, we discuss an experimental extraction of the quantum capacitance from 
experimental data. 

PACS numbers: 73.21.-b, 73.22.-f, 73.20.-r, 81.05.Uw 
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I. INTRODUCTION 

Graphene, a two-dimensional (2D) honeycomb struc- 
ture of carbon atoms, has attracted a lot of interest since 
its isolation in 2004^. It demonstrates unique proper- 
ties which originate from the Dirac-type spectrum of low- 
energy quasiparticles. Nowadays, graphene is considered 
to be a viable alternative to Si for the channel of field- 
effect transistors (FETs) 2 -. One of the main characteris- 
tics of such devices is a capacitance formed between the 
channel and the gate. The capacitance is important for 
understanding fundamental electronic properties of the 
material such as the density of states (DOS) as well as 
device performance including the I — V characteristics 
and the device operation frequency. 

In a classical regime, the capacitance describes the ca- 
pability of an object to store electrical charges and is 
completely determined by the object's geometry and a 
dielectric constant of the medium. If the object's size 
shrinks to a nanometer scale, quantum effects have to be 
taken in account. One of manifestations of these effects 
is a finite DOS which originates from the Pauli exclu- 
sion principle. Low-dimensional systems, having a small 
DOS, are not able to accumulate enough charge to com- 
pletely screen the external field. In order to describe 
the effect of the electric field penetration through a two- 
dimensional electronic gas (2DEG) Luryi introduced a 
concept of a quantum capacitance 3 . 

Recently, the quantum capacitance of a bulk graphene 
layer deposited on a gated Si02 insulated surface 
has been investigated by means of scanning probe 
microscopy^. To the best of our knowledge, no stud- 
ies of the gate capacitance of the graphene nanoribbons 



(GNRs) have been reported yet. However, such studies 
are already technologically feasible. Indeed, during last 
years the great progress has been achieved in fabrication 
and patterning of the GNRs^jSiZiS as well as in controlling 
the morphology, geometry and stability of the graphene 
edges&i£. On the other hand, the quantum and classical 
capacitance of related structures, - carbone nanotubes, 
has been measured and analyzed by a number of groups 
during the last year a n i 12 ' 13 i 14 . The later studies have re- 
vealed a number of interesting properties of the system at 
hand including the structure of the DOS and signatures 
of the electron interaction and correlation. 

In order to provide physical insight into the gate elec- 
trostatics and capacitance of the GNRs, it is important 
to develop intuitive analytical models capturing the es- 
sential physics of the device at hand. Such models are 
also imperative in experimental measurements because 
the quantum capacitance is not directly accessible in the 
experiments and can only be indirectly extracted from 
the measured total capacitance. In the present paper 
we develop a basic analytical theory for the gate elec- 
trostatics and the classical and quantum capacitance of 
the GNRs. We complement this analytical theory by 
exact self-consistent numerical calculations based on the 
tight-binding p-orbital Hamiltonian within the Hartree 
approximation. We demonstrate that the analytical the- 
ory is in a good qualitative (and in some aspects quanti- 
tative) agreement with the exact calculations. There are 
however some important discrepancies. In order to un- 
derstand the origin of these discrepancies we investigate 
the self-consistent electronic structure and charge den- 
sity distribution in the nanoribbons and relate the above 
discrepancy to the inability of the simple electrostatic 
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model to capture the classical gate electrostatics of the 
GNRs. In turn, the failure of the classical electrostatics is 
traced to the quantum mechanical effects leading to the 
significant modification of the self-consistent charge dis- 
tribution in comparison to the non-interacting electron 
description. 

It should be noted that particular aspects of the self- 
consistent gate electrostatics and the electron structure 
of the GNRsi 5 - and the numerical^ and analytical^ stud- 
ies of quantum capacitance of the GNRs have been re- 
ported in the literature. In particular, the quantum 
capacitance of the GNRs as a function of the Fermi 
energy, Cq = Cq(Ef), has been studied in Ref. 
Experimentally, however, the dependence of Cq on the 
Fermi energy is not accessible, and we stress that the fo- 
cus of our analytical and numerical analysis is the gate 
capacitance, C = C(V g ), - the characteristic that is mea- 
sured experimentally (V g being the gate voltage). 

The paper is organized as follows. In Sec. II we for- 
mulate the basics of our model of the gated GNRs. The 
analytical treatment of the gate electrostatics and the 
quantum and classical gate capacitance of GNRs is given 
in Sec. III. The results of the self-consistent numeri- 
cal calculations and a comparison between the analytical 
and numerical calculations are presented and discussed 
in Sec. IV. Section V summarizes the main conclusions. 



II. MODEL 

In experiments graphene samples are separated from 
the gate by a relatively thick insulating substrate (of a 
typical width of at least 300 nm) in order to enable visual 
identification of the graphene sheet. This simplest exper- 
imental setup (with a single back gate) is not particularly 
suitable for measurements of the quantum capacitance 
because in this case (as we will demonstrate below) the 
total capacitance is completely dominated by the clas- 
sical contribution and can hardly be extracted from the 
measured total capacitance. In order to distinguish the 
quantum contribution the gate should be placed much 
closer to the ribbon such that the classical capacitance 
Cc becomes comparable to the quantum one. In our 
study we therefore consider an embedded top-gate ge- 
ometry shown in Fig. [T] where a graphene ribbon of the 
width w is placed on a thick dielectric layer and covered 
by the second much thinner layer of the width d separat- 
ing it from the top gate with the applied gate voltage V g . 
This top-gate geometry was used by Hani et alJ^ and Na- 
tori et ali^ for measurements of the quantum capacitance 
of carbone nanotubes. We also assume that the GNR is 
connected to the source and drain electrods playing a 
role of ideal reservoirs supplying electrons to the ribbon. 
The experimental setup might also include the back gate 
which can, independently of the top gate, adjust a posi- 
tion of the Dirac point in the graphene nanoribbon and 
hence change its electron density. However we assume 
that the back gate is situated much further apart from 



(a) 



top gat 



graphene nanoribbon 
insulating substrate 




12 34 



FIG. 1: (Color online) (a) A schematic diagram illustrating a 
top-gate geometry where an infinitely long graphene nanorib- 
bon of the width w is embedded in a gate insulator with the 
relative dielectric constant e r ; d is the distance between the 
ribbon and the top gate. It is assumed that the graphene 
nanoribbon is connected to the source and drain reservoirs 
supplying electrons to the ribbon, (b) An armchair graphene 
ribbon of the width N. All the results presented in this paper 
correspond to N — 98 (w = 12 nm.) 



the ribbon in comparison to the top gate, and therefore 
for the sake of simplicity we in our model disregard a 
possible capacitative coupling between the ribbon and 
the back gate. 

In this paper we considered two representative struc- 
tures, (I) Hf02 insulating layer with d = 30 nm and e r = 
47 m 1 iE 1 20 i an d (n) Si0 2 insulating layer with d = 300 
nm and e r = 3.9. For the first structure Cc is compa- 
rable to Cq, whereas for the second one Cc <C Cq. We 
limit our calculations to the case of the armchair GNRs, 
whereas we expect that the main results and conclusions 
presented in this paper can be extended to the case of the 
zigzag GNRs. (Note that we do not focus on any specific 
peculiarities of the DOS near the Dirac point like surface 
states in the case of zigzag ribbons) . The spin effects and 
the effect of disorder are outside the scope of our paper 
and are deferred to future studies. All results correspond 
to the metallic armchair GNRs with the width w = 12 
nm (N = 98). We also made calculations for a semicon- 
ductor armchair GNR as well as for wider ribbons (w =50 
nm) and all the results show the same features. 

The system presented on Fig. [T] is described by the 
standard p-orbital tight-binding Hamiltonia n 21 ' 22 



E** 



(1) 



where i r , r +A = 2.5 eV is a nearest-neighbor hopping in- 
tegral; Vh (r) is a Hartree potential at the site r which re- 
sults from the Coulomb interaction between extra charges 
q(r) in the system (including the mirror charges) 15 ' 23 ' , 



Vh(t) 
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47r£n£ r , 

r 



v /|r-r'| 2 +4rf 2 / 

The summation in Eq. (|2|) can be split into two parts 
corresponding to the A and B sublattices of graphene 
(see Fig. [1] (b)). Changing summation to integration in 
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the z-direction we obtain, 



N 

E 



(3) 



x In • 



V ij')~ 

where qA(B)W) is the charge on the carbon atom which 
is located on the j' line and corresponds to the A{B) 
sublattice. 

We solve Eq. fTj numerically to find the Green's func- 
tion using the technique described by Xu et alM. This 
technique greatly facilitates computation speed since it 
does not require self-consistent calculation of the surface 
Green's function. The Green's function in the real-space 
representation, G(r, r), provides an information about 
the local density of states (LDOS) at site r, 



p(r,E) = -— 9f[(G(r,r))], 

7TO 



(4) 



where factor 2 indicates a spin degeneracy and S is the 
area corresponding to one carbon atom. The LDOS can 
be used to calculate the local electron density at the site 



r. 



n(r, Ef 



e f 



dEp(r, E)f FD (E-E F ) 



(5) 



eV c 



where Ef — eV g is Fermi energy and /fd is the Fermi- 
Dirac distribution function. (All the calculations rc- 
poreted in this paper correspond to the temperature 
T = K). The position of the charge neutrality point 
eVc at a given gate voltage V g is determined numerically 
from the calculated dispersion relation. For example, for 
the armchair GNRs, the position of the charge neutral- 
ity point eVc corresponds to the energy which gives the 
minimum number of propagating states with the smallest 
absolute value of the wave vector. Note, that in order to 
achieve a fast convergence, the itegration in Eq. ([5]) is 
performed in a complex plane, since on the real E-&xis 
p(r, E) is a rapidly varying function of the energy (see 
Refs^^ for details) . 

Since the Hartree potential Vh © depends on the elec- 
tron density n(r) which is a solution of the Schrodinger 
equation with the Hamiltonian |T]) , these equations need 
to be solved iteratively. The iteration process is executed 
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until the convergence criterion is met, ym +V m 

where and V™ t are the input and output average 
values of the Hartree potential on the m-th iteration. In 
order to accelerate convergence we used the Broyden's 
second method- 6 , which allows us to reduce the num- 
ber of iterations to ~ 8 — 10 in comparison to ~ 40 — 50 
iterations needed with the "simple mixing" method. 

Having calculated the electron density and the posi- 
tion of the Dirac point numerically, we are in position to 
find the total, quantum and classical capacitances as a 
function of the gate voltage, see Sec. IV for details. The 
analytical approach to the quantum and classical capac- 
itance of the GNRs is described in the next section. 



0.6 

0.5 

0.4 

> 0.3 
03 



LU 



0.2 
0.1 
0.0 
-0 1 

0.3 
0.2- 

> 01 " 

A o.o 



LU 



-0.1 
-0.2 
-0.3 



Va 



Mo 



-0.1 -0.05 0.0 0.05 0.1 

k(1/a) 

FIG. 2: (Color online) A diagram illustrating the change in 
the band structure of the graphene nanoribbon and a shift 
of the chemical potential fi under the application of the gate 
voltage V g (see text for details). The displayed diagram cor- 



responds to the Hf02 structure with d - 
V. 



30 nm and V„ = 0.5 



III. CLASSICAL AND QUANTUM 
CAPACITANCE OF GRAPHENE 
NANORIBBONS: AN ANALYTICAL MODEL 

The electronic structure of planar gated graphene 
sheets was studied by Fernandez-Rossier et alJ^. In this 
section we follow the approach outlined by Fernandez- 
Rossier et al^ and provide an analytical description of 
the one-dimensional charge density and quantum and 
classical capacitance of the graphene nanoribbons. 

Application of the gate voltage Vg to a metallic gate 
induces extra carriers with the density n to the nanorib- 
bon as well as extra carriers of the opposite polarity to 
the gate itself. (Note that the gate and the ribbon repre- 
sent together a charge neutral system). The application 
of the gate voltage shifts the chemical potential \i of the 
ribbon from the charge neutrality point 



eV g = n - no, 



(6) 



It is convenient to represent this shift as a sum of two 
terms 



A* - Mo = eVc + eV Q , 



(7) 



where eVc describes the position of the charge neutrality 
point at the applied voltage V g , and eVQ describes the 
change in the chemical potential due to the filling of the 
quantum mechanical energy bands, see Fig. [2] for illus- 
tration. For a classical conductor the density of states 
is infinite and thus Vq — 0. This provides a natural in- 
terpretation of Vc as a classical electrostatic potential, 
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whereas the potential Vq has a quantum mechanical ori- 
gin reflecting the structure of the quantum mechanical 
density of states. Using a relation 



Vg=V Q + V C 



(8) 



which follows from Eqs. ([6]) and ([7]) and using a definition 



of a capacitance C = 7^7-, we obtain^ 



av 

— 1 /~1 — 1 

tot — U C 



(9) 



g^r, and the classi- 



where the total capacitance Ct t 
cal and the quantum mechanical capacitances are respec- 
tively C C = ff§ and C Q = §f§. 

In order to find the total capacitance for a given 
gate voltage we have to calculate the electron density 
n. The later at the zero temperature is given by n = 
f£ p{E)dE - fj£ p {E)dE 1 where p (E) is the density 
of states for a charge neutral ribbon (fi — (iq), p(E) is 
the DOS at fx = u + eV g , and the integration starts from 
the bottom of the energy band Ej,. (In the following we 
will refer to the case V g = as to an uncharged ribbon, 
and to the case of V g ^ as a charged ribbon). Neglect- 
ing changes in the DOS of the charged ribbon under the 
applied gate voltage in comparison to the uncharged one, 



p(E) =p (E-eV c ), 
the density of the extra carriers reads^ 



Vo+eV Q 



Po(E)dE. 



(10) 



(11) 



(In the next Section we will demonstrate that this ap- 
proximation holds extremely well despite of some modi- 
fications of the band structure for higher gate voltages). 

The DOS of graphene armchair nanoribbons can be 
written in the form (see Appendix I) 



Po(E) 



E 



\E\ 



7rV3ta ^ y/E 2 - El 



\E\-\E n \), (12) 



where n = 0, ±1, ±2, . . . , t w 2.5 eV is the first-neighbor 
hopping integral, a = 0.246 nm is the graphene lattice 
constant (note that = VfH with vf being the Fermi 
velocity), and E n are the subband threshold energies 
whose analytical expressions are provided by Onipkc 27 
(see Appendix I for the explicit expressions for E n ). Sub- 
stituting this expression of the DOS into Eq. (TlTj) , we 
obtain for the one-dimensional (ID) electron density of 
the nanoribbon, 



E 



eV Q y-E*9(\eV Q \-\E n \), 



(13) 

where n — 0, ±1, ±2, .... Here and hereafter without loss 
of generality we set liq = 0. Using the definition of the 
quantum capacitance we get 

edn 



fl7 



dV 



e Po(eV Q ). 



(14) 



In order to calculate Cq from the above equation we 
have to know the position of the chemical potential with 
respect to the charge neutrality point of the charged rib- 
bon, eVg. This can be done from Eq. ([5]), where the 
classical electrostatic potential can be easily calculated 
using the standard method of images, 



V c = 



TVS 



2d arctan 



4d 




(15) 



(the definitions of w and d are given in Fig. [T]). In the 
derivation of this expression we assumed that the surface 
charge density of the graphene nanoribbon a = n/w is 
constant (as expected for a classical capacitor). We will 
discuss the validity of this assumption in the next Sec- 
tion. The classical capacitance (per unit length) of the 
graphene nanoribbon follows from Eq. (1151) 



edn 
dVc 



new 



2d arctan 




Note that in the limit of a narrow ribbon, w <C 
above expression simplifies to Cc = 27re/ln— . 

To summarize, Eqs. ©, ©, ([13 )) - ([TB" )) provide the 
analytical expressions for the ID electron density and the 
total, quantum and classical capacitances of the graphene 
nanoribbons. In order to express the density and the 
capacitances as a function of the gate voltage V g (rather 
than Vq which is not accessible experimentally) we first 
choose some value of Vq and calculate n and Cq from 
Eqs. (fT3|) . (fT4"|) . We then use the calculated values in 
Eqs. (fl"5|) , (fTBj) to find corresponding Vc and Cc- Finally, 



relating Vq to V g via Eq. ([5]), we express the density and 
the capacitance as a function of the gate voltage V g 



IV. RESULTS AND DISCUSSION 

Figure [3] shows the analytical and numerical densi- 
ties and capacitances for two representative nanoribbon 
structures introduced in Sec. II. The left column cor- 
responds to the 30 nm Hf02 dielectric structure with 
e r = 47, and the right column corresponds to a conven- 
tional 300 nm Si02 structure with e r = 3.9. The analyti- 
cal results are based on the expressions given by Eq. ©, 
([9]), (fl"3j) - (fl"6]) and show the electron density n, the total, 
quantum and classical capacitances Ctot,CQ,Cc (Figs. 
|3](a), (b), (c), (d) respectively). 

The numerical results are based on the self-consistent 
solution of Eqs. CQ)-(© as described in Sec. II. We first 
calculate the electron density as a function of the gate 
voltage and then differentiate it numerically in order to 
compute the total capacitance, C tot = see Figs. [3] 
(a) and (b) respectively. Figure [3](c) shows the quantum 
capacitance Cq which is calculated from the DOS on the 
basis of Eq. dHJ) (Note that p {E) = T, r p{ r ,E), 
where p(r,E) is the LDOS given by Eq. 0]), and sum- 
mation is performed over one unit cell containing N ce u 
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FIG. 3: (Color online) The analytical and numerical depen- 
dencies on the applied gate voltage of (a) the electron density 
n, (b) the total capacitance Ctot, (c) the quantum capaci- 
tance Cq, (d) the classical capacitance Cc- (c) also shows 
the "extracted" quantum capacitance (see text for details). 
Left and right panels corresponds respectively to Hf02 and 
Si02 structures. 



sites). Having calculated Ctot and Cq we compute the 
classical capacitance from Eq. © as C^ 1 = C t ~* — Cq , 
see Fig. [3] (d). 

The total capacitance C to t for both structures shows 
characteristic features that can be traced to the corre- 
sponding features in the quantum capacitance Cq, (cf. 
Figs. [3] (b) and (c)). Because the quantum capacitance 
is proportional to the DOS (see Eq. (|14p) the peaks in Cq 
signal consecutive population of electron subbands as the 
gate voltage increases. Note that these features in Ctot 
are much less pronounced for the case of a conventional 
SiC>2 structure because its classical capacitance is much 
smaller than the quantum one (Cq/Cc ~ 40), whereas 
for the HTO2 structure this ratio is only ~ 2 (note that 
the quantum and classical capacitances are added in se- 



ries, Eq. ©). 

The comparison of the analytical and numerical cal- 
culations demonstrates that the analytical theory qual- 
itatively reproduces the exact results very well. There 
is however some quantitative discrepancy in the values 
of both quantum and classical capacitances. In particu- 
lar, the analytical classical capacitance differs by 20-35% 
from its exact numerical value, see Fig. [3] (d). For the 
case of the SiC>2 structure the numerical Cc shows a slow 
increase as V g increases. This is in apparent contrast with 
the behaviour of its analytical counterpart which is inde- 
pendent on the applied voltage. As far as the quantum 
capacitance is concerned, a visual inspection of Fig. [3] 
(c) indicates that the analytical and numerical Cq would 
coincide if one stretches the scale of V g for the analytical 
capacitance (alternatively contracts the scale of V g for 
the numerical capacitance). 

In order to understand the differences between the an- 
alytical and the exact results let us critically inspect the 
assumptions that have been made in the derivation of 
the analytical expressions in the previous section. Let us 
start with the classical capacitance Cc given by Eq. (|16|) . 
In its derivation we assumed that the induced charge 
density is homogeneous and the potential of the ribbon 
is constant as expected for a classical conductor. Fig- 
ures ([!]) (b),(c) show respectively the electron density 
distributions and the Hartree potential for the various 
nanoribbon structures for different values of Vq and thus 
for different densities n (Note that the amount of the in- 
duced charge density n is completely determined by the 
value of Vq defining the position of the Fermi energy with 
respect to the charge neutrality point, Eqs. (fTTj) . (|13p . 
For a reference purpose, the values of Vq are indicated 
at the corresponding dispersion relations shown in Fig. [4] 
(a) ). In order to outline the role of the electron-electron 
interaction we show both the self-consistent Hartree and 
the noninteracting one-electron calculations (respectively 
right and left parts of the panels in Fig. [4]). 

The electron density distribution in the GNRs shows 
the pronounced oscillation between neighboring sites. We 
therefore also show the electron density averaged over 
several neighboring sites, Fig. 2] (b). For a small gate 
voltage (HfC>2 structure, V g = 0.1V) when only the first 
subband is filled the averaged electron density distribu- 
tion is almost uniform and there is practically no dif- 
ference between the self-consistent and the one-electron 
approaches. By increasing the gate voltage to V g — 0.5V 
more carriers are induced and the electron density in- 
creases near the edges of the structure due to the elec- 
trostatic repulsion. It is important to stress that it is not 
only the concentration of the induced charge but primar- 
ily the applied gate voltage that determine the charge 
density distribution. Figure [4] (b) shows the electron 
density distribution and the potential profile for the SiC>2 
structure for the gate voltage V g — 5.75V. This gate volt- 
age is chosen such that the value of Vq = 0.25 eV is the 
same as for the HfC>2 structure with V g — 0.5V (see Fig. 
2] (a)), i.e. the induced charge concentrations are simi- 
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FIG. 4: (Color online) A comparison between the non- 
interacting case (left parts of the panels) and the self- 
consistent Hartree model (right parts) for the Hf02 structure 
with V g — 0.1 V (first column), HfC>2 structure with V g = 0.5 
V (second column), and Si02 structure with V g = 5.75 V 
(third column), (a) The band structure of the nanoribbons; 
thin solid lines indicate the positions of the Fermi energy. 

(b) The electron density distribution across the nanoribbon. 

(c) The potential profiles in the Hartree approximation (right 
parts); the numerically calculated value of Vc (left parts). 
In (b) and (c) the black lines indicate the values of the elec- 
tron density and the Hartree potential averaged over several 
neighboring sites. 



lar. However, the density distribution profile for the SiC>2 
structure is strikingly different showing a strong redistri- 
bution of the charges toward the edges when the applied 
voltage is increased. The larger the applied voltage, the 
stronger the redistribution of the electron density. This 
explains the observation that the numerical Cc gradually 
changes when the gate voltage is increased (see Fig. E](d), 
right panel) . We therefore conclude that the assumptions 
appropriate for a classical capacitor (the charge density is 
homogeneous and the potential of the ribbon Vc is con- 
stant) are violated for the graphene nanoribbons which 
leads to the difference between the analytical theory and 
the exact numerical calculations. 

Note that the macroscopic charge accumulation along 
the boundaries of the graphene strip was discussed by Sil- 
vestrov and Efetov^. Their semiclassical approach and 
the exact numerical calculations presented here demon- 
strate that the density distribution and the potential pro- 
file in the GNRs are qualitatively different from those 
in conventional split-gate quantum wires with a smooth 
electrostatic confinement where the potential is rather 
flat and the electron density is constant throughout the 



wire22. At the same time, the electron distribution and 
the potential profile in the GNR are very similar to those 
in the cleaved-edge overgrown quantum wires (CEOQW) . 
Indeed, the potential profile in the CEOQWs also ex- 
hibits triangular-shaped quantum wells in the vicinity 
of the wire boundaries and the electron density is also 
strongly enhanced close to the edges^. This similarity 
simply reflects the fact that both the CEOQWs and the 
GNRs correspond to the case of the hard-wall confine- 
ment at the edges of both structures. Let us now discuss 
the quantum capacitance Cq and the difference between 
the corresponding analytical and numerical results. The 
crucial assumption used in the analytical model is that an 
application of the gate voltage V g simply shifts the bands 
on the amount Vc such that the DOS of the GNR remains 
unchanged relatively to the charge neutrality point for 



any V g , Eq. (jlOj) . Because the amount of the induced 
charges is completely determined by the value of Vq 
(Eq. I13[) . this assumption implies that the dependence 
n = ti(Vq) obtained by both analytical and numerical 
calculations should coincide. In order to verify this as- 
sumption we compare the evolution of the band diagram 
for structures with different classical capacitances for dif- 
ferent gate voltages. Figure[5](a) shows the dispersion re- 
lation for Hf02 and Si02 structures where the gate volt- 
ages V g are chosen such that Vq is the same in both cases 
(V Q = 0.25V for V g m ° 2 = 0.5V and V g si ° 2 = 5.75V). (For 
comparison we also display the dispersion relations for 
non-interacting electrons). Even though the equal val- 
ues of Vq imply the same induced charge density, the 
changes of the dispersion relations are quite different. 
In the considered V g interval the dispersion relation of 
the HTO2 structure remains practically unchanged, and 
therefore the analytical and the numerical dependencies 
for n — ti(Vq) as well as for Cq = Cq(Vq) are almost 
undistinguished (see Fig. [5](b), (c), left panel). Modifica- 
tion of the band structure is much stronger for the case of 
the Si02 structure with smaller Cc ■ This is also reflected 
in the analytical and numerical dependencies n = ti(Vq) 
exhibiting a difference up to 15% in the considered gate 
voltage interval. However, despite of this difference for 
n = ti(Vq) the corresponding difference between the an- 
alytical and numerical results for Cq = Cq(Vq) is prac- 
tically negligible even for the Si02 structure, see Fig. [5] 
(b), (c), right column. 

The reason for the modification of the band structure 
can be understood from the analysis of the potential dis- 
tribution shown in Fig. [5] (c). In order to shift Vq on 
the same value one should apply a higher gate voltage to 
the structure with a smaller classical capacitance. Differ- 
ent gate voltages applied to different structures produce 
the same shift of Vq but give rise to different distribu- 
tions of the electrostatic potential across the nanoribbon. 
The difference between the electrostatic potential at the 
middle of the ribbon and at the edges, AVh, is an or- 
der of magnitude higher for the structure with smaller 
Cc- Since the ribbon width w is the same for both struc- 
tures, this leads to higher effective transverse electric field 
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FIG. 5: (Color online) (a) The dispersion relations for the 
non-interacting and the Hartree electrons (left and right parts 
of the diagrams respectively). For the Hartree case the disper- 
sion relation is calculated for Vq = 0.25V; the corresponding 
values of V g are indicated above the figure, (b) Analytical 
(solid lines) and numerical (dashed lines) electron concentra- 
tion n and (c) quantum capacitance Cq as a function of the 
Fermi energy level with respect to the Dirac point (Vq). The 
dotted lines serve as a guide to compare the analytical and 
numerical electron densities calculated for Vq = 0.25V. The 
insets in (b) show the dependence of Vq on the applied gate 
voltage V g . Left and right panels corresponds respectively to 
Hf02 and Si02 structures. 



E 



AVj 



w/2 



f- for the case of SiC"2 structure, which, in turn, 



modifies the band structure of the GNR. Note that the 
effect of the electric field on the band structure and the 
DOS of graphene nanoribbons was reported before by 
many authors, resulting in e.g. the energy-gap modula- 



tion semiconductor armchair ribbons^ or opening of the 
energy gap for the case of zigzag ribbons 3 ^. 

We demonstrated above that the analytical and the 
numerical dependencies Cq — Cq(Vq) show an excellent 
agreement in the considered gate voltage intervals both 
for Hf02 and Si02 structures (i.e. the utilization of Eq. 
(fTU)l is fully justified). However, the analytical and nu- 
merical dependencies Vq = Vq(V 9 ) show some discrep- 
ancy, see insets to Fig. [5] (b). This, in turn, leads to 
a discrepancy between the analytical and numerical Cq 
as a function of V g as shown in Fig. [3] (c). This discrep- 
ancy is manifest itself in the difference of the gate voltage 
scale and not in the difference of the magnitudes of Cq . 
Using Eqs. j8|) and ((9]), the change of Vq can be easily 

related co the change in V g , = %^ = f 1 + 

Because for a given Vq the analytical and numerical Cq 
are practically the same, the difference between the ana- 
lytical and numerical dependencies ^p?- is primarily due 
to the difference of the corresponding classical capaci- 
tances Cc- We therefore conclude that the discrepancy 
between the analytical and numerical Cq as a function 
of V g is related to the difference in the corresponding 
classical capacitances. 

Let us now discuss experimental determination of the 
quantum capacitance Cq . In our both analytical and nu- 
merical approaches we are in position to calculate Cq 
directly. In contrast, Cq is not directly accessible in ex- 
periments. For example, for the case of carbon nanotubcs 
it is the total capacitance Ctot that is measured experi- 
mentally. The quantum capacitance Cq is then extracted 



Q — '-'tot 



a 



C ' 



from Ctot according to Eq. ©, 
where Cq is a corresponding analytical expression for 
the classical capacitance of a nanotube (i.e. a classical 
capacitance between a metallic cylinder and an infinite 
plane) 11 ' 13 . Our calculations presented above demon- 
strate that for the graphene nanoribbons the numerical 
Cq differs from its classical analytical expression given 
by Eq. (fTl)]) . Therefore, a question arises, to what extent 
one can rely on the above procedure for the extraction of 
the quantum capacitance? 

In order to answer this question let us assume that 
our numerically calculated Ctot (shown in Fig. [3] (b)) 
corresponds to the experimental data. We then assume 
that the classical capacitance of the GNR is given by 
the analytical expression (|16|) describing a capacitance 
between a metallic strip and an infinite plane. Finally, 
the quantum capacitance (which we will call "extracted" , 
^extracted -j j g obtained by subtracting Cc from the "mea- 
sured" Ctot according to Eq. © . The "extracted" quan- 
tum capacitance is shown in Fig. [3J (c). For the case of 
the HTO2 structure (where Cc ^ Cq), the behavior of 
(^extracted is qualitatively similar to that of the numerical 
Cq, even though their values differ significantly. How- 
ever, for the case of Si02 structure (where Cc -C Cq) 
the "extracted" quantum capacitance Cg xtractcd does not 
reproduce the numerical Cq even qualitatively with all 
the features related to the quantum mechanical DOS be- 
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ing completely lost. This is simply related to the fact 
that due to the series addition of the capacitances, for 
the case of Cc <§C Cq the total capacitance is completely 
dominated by the classical one, and the features in Cq 
can be reproduced only when Cc becomes comparable to 
Cq. It is interesting to note that the difference between 
the experimentally extracted quantum capacitance and 
its expected value was detected for the case of the car- 
bone nanotubes and was attributed to the strong elec- 
tron correlation and signatures of the Luttinger liquid 
behaviour—. Our calculations indicate that the origin 
of such deviations can have a rather simple explanation 
related to the inability of a standard electrostatics to re- 
produce quantitatively the classical capacitance of the 
structure at hand. 



V. CONCLUSIONS 



finement at the edges of the structure. 

Finally, we discuss experimental determination of the 
quantum capacitance Cq. We demonstrate that the ex- 
tracted Cq might significantly deviate from its actual 
value given by the density of states of the GNRs. This 
deviation is related to the inability of the standard elec- 
trostatics to reproduce quantitatively the classical capac- 
itance of the structure at hand. 
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In the present paper we develop an analytical theory 
for gate electrostatics and classical and quantum capac- 
itance of graphene nanoribbons (GNRs). We compare 
the analytical theory with the exact self-consistent nu- 
merical calculations based on the tight-binding p-orbital 
Hamiltonian within the Hartree approximation. 

We find that the analytical theory is in a good qualita- 
tive (and in some aspects quantitative) agreement with 
the exact calculations. There are however some impor- 
tant discrepancies. In order to understand the origin 
of these discrepancies we investigate the self-consistent 
electronic structure and the charge density distribution 
in the GNRs obtained from the exact numerical calcula- 
tions. We demonstrate that the assumptions appropriate 
for a classical capacitor (the charge density is homoge- 
neous and the potential of the conductor is constant) 
are violated for the graphene nanoribbons which loads 
to the difference between the analytical theory and the 
numerical calculations. In turn, the failure of the clas- 
sical electrostatics is traced to the quantum mechanical 
effects leading to the significant modification of the self- 
consistent charge distribution in comparison to the non- 
interacting electron pictures. We also show that as a re- 
sult of electron-electron interaction the band structure of 
the GNRs modifies as the applied gate voltage increases. 

Our exact numerical calculations show that the density 
distribution and the potential profile in the GNRs are 
qualitatively different from those in conventional split- 
gate quantum wires with a smooth electrostatic confine- 
ment where the potential is rather flat and the electron 
density is constant throughout the wire. At the same 
time, the electron distribution and the potential profile 
in the GNR are very similar to those in the cleaved- 
edge overgrown quantum wires (CEOQW) exhibiting 
triangular-shaped quantum wells in the vicinity of the 
wire boundaries accompanied by the corresponding en- 
hancement of the electron density close to the edges. This 
similarity reflects the fact that both the CEOQWs and 
the GNRs correspond to the case of the hard-wall con- 



APPENDIX A: DENSITY OF STATES OF 
GRAPHENE NANORIBBONS 

The density of states of a quantum wire (including a 
factor 2 for the spin degeneracy) reads, 



, k± n ) 



(Al) 



where fcy (E) and k± n denote the longitudinal (continu- 
ous) and the transverse (quantized) components of the 
wave vector, respectively. The summation in (| Al|) in- 
cludes all transverse modes which energy E n < E. The 
dispersion relation for the nanoribbon of the width N in 
the low-energy limit close to the Dirac point is given by 
Onipko^i, 



U<T 2 



(A2) 



where a — A, Z corresponds to the armchair (A) and 
zigzag (Z) GNRs. The transverse wave vector k^ n 
is given by different expressions depending on whether 
2 ^ +1 ^ * s m t e g er (metallic nanoribbon) or not (semicon- 
ducting nanoribbon), 



k A 




a, metallic 



rrj (n — i) a, semiconducting 



0,±1,±2, . . .. For the zigzag structures, the trans- 



verse wave vector kj_ n has a form: 
z _ 7r(n+|) 



2V3N 



(A4) 



n = 0, 1, 2, The parameter kj_ n describes the shift of 

n-th dispersion branch minima respect to the Brillouin 
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zone centre, 



h a — 



for a = A 
^kll for a = Z 



(A5) 



are the subband threshold energies. The electron density 
at zero temperature is obtained by the integration of the 
DOS from the charge neutrality point /iq = to the Fermi 
energy, n = J^ F pdE, 



Using Eq. (|J2|) in Eq. (|A"l"j) . we obtain for the DOS of 

the armchair (zigzag) ribbon 



\E\ 



6(\E\-\E n \), (A6) 



n°(E F ) 



£ 



El 



E°*B{\E F \-\E%\). (A7) 



where n — 0, ±1, ±2, . . . for the armchair GNRs and n — 



0,1,2,... for the zigzag GNRs, and and E° 
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